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Abstract 



Properties of magnetic films are studied in the framework of Ising models. In partic- 
ular, we discuss critical phenomena of ferromagnetic Ising films with straight lines 
of magnetic adatoms and straight steps on the surface as well as phase diagrams 
of the axial next-nearest neighbour Ising (ANNNI) model for thin films exhibiting 
various spatially modulated phases. 
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Magnetism in thin films has attracted much interest, both experimentally 
and theoretically [1,2]. In this paper, we shall deal, in the framework of Ising 
models, with critical properties of films having different types of surface defects 
as well as with phase diagrams of thin films displaying spatially modulated 
magnetic structures, extending our recent work [3,4]. 

To analyse the role of surface imperfections at criticality, we consider the 
ferromagnetic nearest-neighbour Ising model on the simple cubic lattice for 
films of L layers, using Monte Carlo and the density-matrix renormalization 
group techniques. To be specific, we study films with a straight line of magnetic 
adatoms or a straight step of monoatomic height. Both types of defects, lines 
and steps, are assumed to run across the entire surface layer through its center 
and along one of the crystallographic axes. The magnetic couplings at the 
surface, J s , may be identical or different from those in the bulk, J b . The 
local interactions at the defect may be modified, e.g., the couplings between 
neighbouring spins in the additional line [3]. 

Let us first discuss results for films with lines of magnetic adatoms (our work 
has been motivated by related experiments [5] , albeit there the lines of adatoms 
are not coupled magnetically to the substrate).- The temperature, at which 
the transition between the ferromagnetic and the paramagnetic phases occurs, 
rises with the film thickness, T c = T C (L), as shown in Fig.l. Of course, the 
magnetization in the additional line, mi, is usually different from the surface 
magnetization, see also Fig. 1, because of the reduced coordination number 
and the, possibly, modified local couplings. 

To elucidate critical properties at the line as compared to those of the perfect 
surface and of the bulk, one may monitor the effective exponent /3 e //- That 
quantity follows from the ansatz for the magnetization m oc t^ eff , where t 
is the reduced temperature \T C (L) — T\/T C (L). For t — > 0, (3 e ff becomes 
the true critical exponent f3 - In accordance with general considerations on 
critical phenomena at defect lines in systems with two-dimensional critical 
fluctuations [6], we find the critical exponent f3i of the magnetization in the 
line of adatoms to be non-universal. In films of finite thickness L, it varies from 
to 1/2, decreasing when strengthening the local couplings at the line, and 
depending only rather weakly on L [3]. In contrast, the critical exponent of the 
magnetization of a perfect surface is 1/8 for all finite values of L. It remains 
to be 1/8 even in the thermodynamic limit, L — > oo, when the surface orders 
at a higher temperature than the bulk ('surface transition') due to sufficiently 
strong interactions in the surface, J s . When the surface and the bulk order 
at the same temperature ('ordinary transition'), the critical exponent of the 
surface magnetization is, in the limit L — > oo, about 0.80; this quite large 
value leads, for finite L, to a pronounced increase of the effective exponent 
away from criticality, both for the magnetization of the surface and in the line 
[3]. That behaviour should be taken into account in interpreting properly, e.g., 
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Fig. 2. Simulated specific heat C of the Ising model with equal couplings and with 
a step dividing the system into two halfs of one and two layers; the bottom layer 
consists of 120 x 120 spins. 
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Fig. 3. Phase diagram of the ANNNI model with 9 layers, in mean-field approxi- 
mation; r = ksT/Jo. Phases are denoted in the usual manner, Refs. 4 and 8. 

experimental data. 

Another type of surface defect is the straight step of monoatomic height, divid- 
ing the film in two halfs of thickness L and L + 1 . Accordingly, we are dealing 
with a composite system displaying two distinct phase transitions at T C (L) 
and T C (L + 1) [3], as signalled, e.g., by the two-peak structure of the specific 
heat, see Fig. 2 (attention may be drawn to related recent experimental ob- 
servations [7]). Because T C (L + 1) > T C (L), the magnetization at the step edge 
is expected to vanish on approach to T C (L + 1) similarly to the magnetization 
at the surface of a two-dimensional Ising model. Thence, the critical exponent 
ist expected to be 1/2 [6], in accordance with our simulational findings [3]. 

Secondly, we study the influence of the film thickness on spatially modulated 
magnetic phases in the framework of the ANNNI model [8] on the simple cu- 
bic lattice with L layers [4]. Neighbouring spins in each of the L layers are 
coupled ferro magnetically, Jo > 0. Perpendicular to the layers ferromagnetic 
interactions, J\ > 0, between neighbouring spins in adjacent layers, compete 
with antiferro magnetic couplings, Ji < 0, between axial next-nearest neigh- 
bour spins. Our work extends a recent analysis on the influence of the film 
surface on critical properties of the Lifshitz point [9] . 

For each film thickness L, we determine, using mean-field theory, Monte Carlo 
simulations and low-temperature series, a distinct phase diagram (setting 
Jo = Ji), m the (k = —J2/J1, ksT/Jo) plane, consisting of various modu- 
lated phases. The cases L — 5, 7, and 8 have been depicted before [4]. In 
Figs. 3 and 4, we show the mean-field phase diagrams for L = 9 and 10. 
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Fig. 4. Phase diagram of the ANNNI model with 10 layers, in mean-field approxi- 
mation. 

Compared to the ANNNI model in the thermodynamic limit, L — ► oo, main 
differences due to the finite film thickness are: (i) The ordered phases occur- 
ring directly below the transition line to the paramagnetic phase are either 
symmetric or antisymmetric about the center of the film. In fact, neighbouring 
phases have different parity, with their number of extrema in the layer mag- 
netization m(z), z = 1, 2, ...L, increasing by one when increasing k. (ii) For L 
odd, the center layer of the antisymmetric variants of those phases is disor- 
dered ('partially disordered phases'), (iii) For L odd, a transition line arises 
from (k — 1,T — 0), in addition to the multiphase point at (k = 1/2, T = 0) 
being present at all L. The phases springing from the multiphase point may 
include structures (having one 4-band or one 5-band of layers with equally 
oriented spins) being not stable in the thermodynamic limit. 

Of course, it is of interest to study the robustness of our findings against modi- 
fications of the model parameters. In any event, results may provide some guid- 
ance to explain possible experimental realizations for thin films with magnetic 
superstructures (here, attention is drawn to recent local-spin-density calcu- 
lations of thin Fe films on Co substrates showing distinct and complex phase 
diagram for various film thicknesses [10]). 
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